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Abstract 

Dijet cross sections in neutral current deep inelastic ep scattering have been mea- 
sured in the range 10 < < 10"^ GeV^ with the ZEUS detector at HERA using an 
integrated luminosity of 38.4 pb^^. The cross sections, measured in the Breit frame 
using the jet algorithm, are compared with next-to-leading-order perturbative 
QCD calculations using proton parton distribution functions. The uncertainties of 
the QCD calculations have been studied. The predictions are in reasonable agree- 
ment with the measured cross sections over the entire kinematic range. 
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1 Introduction 



Dijet production in deep inelastic ep scattering (DIS) provides a test of perturbative 
Quantum Chromodynamics (pQCD) and is sensitive to the structure of the proton. A 
comparison of measurements of dijet cross sections and pQCD predictions tests the vahdity 
of the concept of factorisation [|l|] into a hard partonic cross section and universal parton 
distribution functions (PDFs). The large centre-of-mass energy of the HERA ep collider 
(y^ 300 GeV) permits measurements of cross sections 0, for jets of high transverse 
energy, covering over three orders of magnitude in both the photon virtuality, Q^, and 
the fraction of the proton's momentum carried by the struck parton, x. 

In leading-order (LO) QCD, two processes contribute to dijet production in DIS: the QCD- 
Compton (QCDC) and the boson-gluon fusion (BGF) processes, shown in Figs, ^(a) and 
(b), respectively. The QCDC and the BGF cross sections are calculated in pQCD by 
convoluting the matrix element for the hard process, which depends upon the value of the 
QCD strong coupling constant, as, with the PDFs of the proton. In the high-Q^ region, 
the QCDC process is dominant and the quark PDFs are well constrained by inclusive DIS 
data. Hence the dijet cross-section measurements allow tests of pQCD and a measurement 
of as- The extraction of ag from the high-Q^ dijet data and its evolution with are the 
subjects of a separate publication The BGF process is the dominant contribution to 
dijet production at < 500 GeV^. Therefore, measurements of the dijet cross section 
at low are sensitive to the gluon momentum distribution in the proton, xg{x), at low x. 
The dijet cross section can be compared to the predictions of NLO QCD calculated using 
various parametrisations of the proton PDFs. These comparisons complement analyses 
that determine xg{x) from the scaling violations of the structure functions [^, |^. 

In this paper, the value of as{Mz) is fixed and the inclusive dijet data (> 2 jets) are used to 
test pQCD and the universality of the proton PDFs, in particular the gluon distribution. 
The precision of the next-to-leading-order (NLO) QCD calculations plays an important 
role in this analysis and is also studied here in detail. The high-statistics dijet sample, 
together with improved NLO QCD calculations and a better understanding of jet-finding 
algorithms, permit higher-precision tests over a wider range of Q^, 10 < < 10^ GeV^, 
than in previous publications 0. 



2 Theoretical framework 



Within the framework of pQCD in DIS, the dijet production cross section, da, can be 
written as a convolution of the proton PDFs, fa, with the partonic hard cross section, 

da = ^ dx fa{x, fi%) daa{x, as{^lR),^i\, /i|) ■ (1 + 5had)- (1) 

The partonic cross section describes the short- distance structure of the interaction and is 
calculable as a power-series expansion in the strong coupling constant. As can be seen in 
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Eq. (IID, the cross section depends on the renormahsation scale, hr, since the calculations 
are not carried out to all orders. The PDFs contain the description of the long-distance 
structure of the incoming proton. The evolution of the PDFs with the factorisation scale, 
fip, at which they are determined follows the DGLAP equations The hadronisa- 
tion correction, 5had, can be estimated using Monte Carlo models for fragmentation (see 
Sections || and [^21) . 

The predictions of QCD were calculated at NLO in using the programs MEPJET 
and DISENT 0. These programs yield parton-level cross sections and allow for an 
arbitrary jet-definition scheme with user-defined cross-section cuts. Unless otherwise 
stated, all the pQCD predictions presented in this publication were calculated using the 
CTEQ4M [T^ proton PDFs. The value, as{Mz) = 0.116, and the formula for the running 



of as used in MEPJET and DISENT were the same as those used in the proton PDF de- 
termination. According to calculations using DISENT, the percentage of gluon-initiated 
dijet events in the total dijet cross section increases from ~ 10% to ~ 70% as decreases 
from 10^ to 10 GeV^. 

In all NLO QCD calculations, fip was set to Q^. The dependence of the pQCD predictions 
on fip and yu|. is discussed in Section \lA\ . In DIS dijet events, there are two "natural" 
variables to define the of the event and the square of the transverse energy of the 
event or jets. In principle, any multiple of these variables or any combination of them can 
be used. The interplay between these two variables has been investigated by measuring 
the dependence of the dijet cross section on the ratio E^/Q'^, where Et is the average 
transverse energy of the two jets with highest transverse energy in the event. 



3 Experimental setup 

The dijet data sample presented here was collected with the ZEUS detector during 1996 
and 1997 and corresponds to an integrated luminosity of 38.4 ± 0.6 pb^^. During this 
period, HERA operated with protons of energy Ep = 820 GeV and positrons of energy 
E, = 27.5 GeV. 



The ZEUS detector is described in detail elsewhere |TT|, |T2|. The main components used 
in the present analysis are the uranium-scintillator sampling calorimeter (CAL) |T3|, and 



the central tracking chamber (CTD) positioned in a 1.43 T solenoidal magnetic field. 
The smallest subdivision of the CAL is called a cell. The CAL relative energy resolutions, 
as measured in test beams, are 18%/v/^(GeV) and 35%/^E{GeV) for electrons and 
hadrons, respectively. The interaction vertex is measured using the CTD with a typical 
resolution along (transverse to) the beam direction^ of 0.4 (0.1) cm. 



The luminosity was measured using the Bethe-Heitler reaction e~^p e^'^p |T5|. The 
resulting small-angle energetic photons were measured by the luminosity monitor, a lead- 
scintillator calorimeter placed in the HERA tunnel aX Z = —107 m. 



^The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity is defined 
as 77 = — In (tan ^), where the polar angle, 9, is measured with respect to the proton beam direction. 
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4 Reconstruction of kinematic variables 



The reaction 

e+{k)+p{P) ^e+{k')+X 

at fixed squared centre-of-mass energy, s = {k + P)^, can be fully specified in terms of 
= —q^ = —{k — k'Y and Bjorken x = Q'^/{2P ■ q). The fraction of the positron's 
energy transferred to the proton in its rest frame is y = Q'^/ (sx). 

For processes where two or more jets are produced in the final state, the observable ^ is 
defined as 

X 1+ 



where the dijet mass, Mjj = ■\/2£^ii?2(l — cos0i2), is calculated from the energies Ei and 
E2 and the opening angle 612 of the two jets of highest transverse energy. The variable 
^ is the fraction of the proton momentum carried by the struck parton in the LO QCD 
processes (see Fig. ||). 

The kinematic variables were reconstructed using a combination of the electron and the 



double-angle (DA) methods [|T6|. The electron method was used except when the angle 
of the hadronic system, ■jh, was less than 90° and the scattered-positron track could be 
well reconstructed by the CTD. In this case, the DA method was used. 

The variable yjB = '~ cos 9i) / {2Ee) , calculated according to the Jacquet-Blondel 

method |T^, where the sum runs over all CAL cells except those belonging to the scattered 



positron, gives a measurement of y with good resolution at low y. 



5 Selection of the dijet event sample 



The events were selected online via a three-level trigger system |T^, |19] using the same 



selection algorithms as in previous dijet publications pO|. Neutral current DIS events 



were selected by requiring that the scattered positron was measured in the CAL pT |. 
Further criteria were applied both to ensure an accurate reconstruction of the kinematic 
variables and to increase the purity of the sample: 



• ii^g > 10 GeV, where E'^ is the scattered positron energy, after correction for energy 
loss in inactive material in front of the CAL, to achieve a high-purity sample of DIS 
events; 

• 38 < E — pz < 65 GeV, where E — pz = Ei{l — cos 9i) and the summation is 
over all CAL cells, to remove background from photoproduction and events with 
large initial-state QED radiation; 
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He < 0.95, where i/e = I — — cos6'g) and 9'^ is the polar angle of the scattered 

positron. Along with the previous requirements, this reduces the photoproduction 
background to a negligible level 

|X| > 14 cm or |y I > 14 cm, where X and Y are the impact positions of the positron 
on the CAL, to avoid the low-acceptance region adjacent to the rear beampipe; 

l-Z^vertexl < 50 cm, to eusure that event quantities can be accurately determined; 

yjB > 0.04, to give sufficient accuracy for DA reconstruction of and x; 

10 < g2 < GeV2. 



After these cuts, jets were reconstructed using the longitudinally invariant kx cluster 
algorithm p3 in the inclusive mode. The jet search was conducted in the Breit frame 



defined by q + 2xP = 0, where q and P are the three- momentum vectors of the exchanged 
boson and the proton, respectively. In the Breit frame, the three-momentum vectors of 
the exchanged boson and the incoming and outgoing partons all lie on the Z axis for a 
quark-parton-model type of event, i.e. the production of a single quark in the final state. 
For a typical QCDC or BGF event, the transverse energy in this frame is non-zero. The 
vector required to boost to the Breit frame was determined using the reconstructed event 
kinematics. 

For each event, the jet search was performed over all the CAL energy deposits, considered 
as massless objects and boosted to the Breit frame, excluding those corresponding to the 
identified scattered positron candidate. The jet transverse energies were corrected for 
energy loss in the inactive material in front of the CAL. Events with two or more jets 
found in the Breit frame were selected by requiring that the two jets with the highest 
transverse energy satisfied the following cuts: 



• E^^^ > 8 GeV and > 5 GeV, where E^^'^ and E^^" the transverse 
energies in the Breit frame; 

• E^^^ > 5 GeV and -E^^^ > 5 GeV, where EJj^^^ and E^^^ are the transverse 
energies of the two jets in the laboratory frame; 

• l'7i'2^^l < 2, where ri]'2^ are the pseudorapidities of the two jets in the laboratory 
frame. 



The jet transverse energy and pseudorapidity cuts in the laboratory frame were imposed in 
order to select well-measured jets within the acceptance of the CAL. After all cuts, 39576 
events with two or more jets, including 3902 events with three or more jets, remained in 
the sample. 
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6 Monte Carlo simulation 



Monte Carlo (MC) simulations were used to correct the data for detector acceptance and 
resolution. Two MC models were used to generate DIS events: ARIADNE 4.10 |^ and 
LEPTO 6.5 



26 



In ARIADNE, the QCD cascade is simulated using the colour-dipole 
model. LEPTO uses the exact matrix elements to generate the hard process and the 
parton-shower model to simulate higher-order processes. Both models use the Lund string- 



fragmentation model ||2^ for hadronisation, as implemented in JETSET 7.4 . To take 



into account first-order electroweak corrections, LEPTO and ARIADNE were interfaced 



with HERACLES 4.5.2 [|| using the DJANG06 2.4 Q program. The CTEQ4D proton 
PDFs were used. To estimate the uncertainty due to hadronisation, events were also 
produced using the HERWIG 5.9 generator in which the fragmentation is simulated 
using a cluster model 



The ZEUS detector response was simulated with a program based on GEANT 3.13 p3 . 
The generated events were passed through the simulated detector, subjected to the same 
trigger requirements as the data and processed by the same reconstruction and offline 
programs. 

Figures |] and |] show normalised uncorrected differential distributions of dijet events 
together with the predictions of the ARIADNE and LEPTO MC programs. As shown 
in Fig. 0(a), (b) and (d), the global event variables are better described by ARIADNE, 
while both ARIADNE and LEPTO give an adequate description of the data in Fig. 0(c). 
For the jet quantities, LEPTO gives a better description of the data in Fig. 0(a), while 
for the three other distributions, ARIADNE is better overall. Therefore, ARIADNE was 
used as the default MC simulation to determine the corrections from the detector to the 
hadron level. 



7 Theoretical uncertainties 



7.1 Uncertainties due to the renormalisation and factorisation 
scales 

A study of the uncertainties introduced in the theoretical predictions from the renormal- 
isation scale has been performed by choosing to be either or where is 
the square of the sum of the transverse energies of all final-state partons, and by varying 
each by a factor of four. 

The ratio of the predicted dijet cross sections, o"(/i|j = Q'^ / A) / a{^\ = Q"^), as a function of 
and calculated using DISENT is shown in Fig. 0(a). The theoretical uncertainty 
depends mainly on and decreases from values of about 50% at low to about 10% at 
high Q^. The ratio of the predicted dijet cross sections, cT(/i|. = £'|./16)/cr(/x|j = E^/i), 
as a function of and E^/A calculated using DISENT is shown in Fig. 0(b). The 
estimated theoretical uncertainty observed in Fig. 0(b) is always less than that shown in 
Fig. 0(a) and has a similar dependence on both and E^/A over the entire kinematic 
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range. Results similar to those shown in Fig. |^ were obtained when calculating the ratios 
cr(/x| = 4Q2)/cr(/i| = Q^) and a{fi% = E^)/a{fi% = Even for Et as high as 

100 GeV, the scale uncertainty of the pQCD calculations is smaller than 10% only for 
greater than 2000 GeV^ 

The uncertainties introduced in the theoretical predictions from the factorisation scale 
were studied by changing the value of yu|. = by a factor of four. The resulting change 
in the cross section was about 5%, significantly less than when /i^ was varied by the same 
factor. 

7.2 Hadronisation uncertainty 

The results of the NLO QCD programs for calculating dijet production are given in terms 
of jets of partons, while the measured cross sections refer to jets of hadrons. In both 
cases, the jets were reconstructed using the longitudinally invariant fcr algorithm in its 
inclusive mode. Therefore, a correction was applied to the predicted parton-level cross 
sections to account for the effects of hadronisation. Correction factors were defined for 
each bin as the ratio of the parton-level to hadron-level cross sections and were calculated 
using the MC simulations. The predictions of NLO QCD were divided by the mean of 
the correction factors calculated using the ARIADNE and LEPTO MC programs. Since 
LEPTO with default parameters does not agree with the data as well as ARIADNE (see 
Fig. 0), the LEPTO sample was reweighted to agree with the distribution of the data. 
After this reweighting, all other kinematic distributions were found to agree better with 
the data. The difference in the LEPTO correction factors before and after reweighting 
was negligible. 

The correction factors as a function of were found to be ~1.25 at low Q^, falling 
to ~1.15 for > 100 GeV^. As a function of log^oC and E^/Q"^, the correction was 
in the range 1.15 — 1.30. The difference between the mean correction factor and the 
individual correction factors from LEPTO and ARIADNE was defined as the uncertainty 
of the theoretical predictions due to uncertainties in the parton-to-hadron corrections. On 
average, this difference is about 5%. It is largest (12%) in the range —2.4 < log^g ^ < ~2.2, 
while, for < 25 GeV^, a difference of about 10% was observed. The correction factors 
obtained from HERWIG, which uses a different hadronisation model than either LEPTO 
or ARIADNE, were also calculated and were found to deviate by typically less than 5% 
from the correction factors obtained with ARIADNE. 



8 Data corrections and systematic uncertainties 

The cross sections for jets of hadrons in bins of log^gQ^^ logio^ 1oSio(-^t/'5^) were 
obtained by applying a bin-by-bin correction to the measured dijet distributions using 
ARIADNE. The corrections take into account the efficiency of the trigger, the selection 
criteria and the purity and efficiency of the jet reconstruction. An additional MC cor- 
rection was applied to the measured cross sections to account for initial- and final-state 
QED-radiation effects. 
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A detailed study of the systematic uncertainties of the measurements was performed p4| , 
|35[] . The uncertainty due to that of the jet-energy scale was studied by selecting a separate 
sample of DIS events with only one jet in the final state in the laboratory frame. In 
this sample, the transverse energy of the scattered positron, which is known to within 



±1% |2T[, balanced the transverse energy of the jet. The ratio between the jet and 
the positron transverse energies was evaluated as a function of the jet pseudorapidity, 
separately for data and MC-generated events. A difference of up to ±3% was observed, 
resulting in a systematic uncertainty of about ±10% for the measured cross sections 
throughout the kinematic range. The jet-energy scale dominates all other uncertainties 
except at > 1000 GeV^, where the statistical plus other systematic uncertainties 
become comparable. 

The other systematic uncertainties originate from the residual uncertainties in the event 
simulation and were estimated by correcting the data using the LEPTO correction factors 
and by varying the cuts in both data and MC simulation by an amount equal to the 
resolution on the relevant quantity. The most important were: 



• use of the MC model - using LEPTO instead of ARIADNE to evaluate the accep- 
tance corrections resulted in an uncertainty of typically ±5%; 

• jet-transverse- energy cuts - E^^^, E^^^, EJj^^^ and EJj^^^ were simultaneously varied 
by the jet-transverse-energy resolution near the cuts, ±17%, to account for differ- 
ences between the data and the MC simulation. This resulted in an uncertainty of 
up to ±12% and typically of ±4%; 

• jet-pseudorapidity cuts - a change of ±0.1 in the //i^^ cuts imposed on the jets in 
the laboratory frame resulted in an uncertainty of up to ±6% and typically of ±2%; 

• E — Pz cuts - the lower E — Pz cut of 38 GeV was varied by ±4 GeV, resulting in 
an uncertainty of less than ±2%. 



All the uncertainties are correlated between bins except for that due to the acceptance 
correction using LEPTO. Nevertheless, they have been added in quadrature when dis- 
played in the figures and tables, with only the uncertainty due to the jet-energy scale 
shown separately. 



9 Results and discussion 



9.1 Asymmetric jet cuts for the dijet cross section 

Figure H shows the measured and predicted inclusive dijet cross sections at the hadron 
level as a function of the threshold on the transverse energy of the leading jet, E'^^^''^"*. 
For this study, the requirement on the jet with the highest transverse energy was relaxed 
and the phase space region was defined by: 10 < < 10^ GeV^, y > 0.04, E'^ > 10 GeV, 
E^^^ > 5 GeV, E^^^ > 5 GeV, E^^b > 5 QeV, E^ab > 5 QeV and \r]^^^\ < 2. 
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The NLO QCD predictions of DISENT with /i^ = exhibit an unphysical behaviour at 
low E^^'^"'^; as the threshold decreases below 6 GeV, the predicted cross section 

decreases, whereas the dijet cross section should increase. This occurs because, within 
NLO QCD near the "symmetric-cut" threshold at which E^^'^^^ = E'^^^'"^"*, the phase 
space for the real emission of a third parton is reduced. This results in an incomplete 



cancellation between the real emission and the virtual-loop corrections |^ and hence the 



cross section falls. As seen in Fig. & the predictions of DISENT with jjj^ = Q approach 
the measured cross section for Ej,^ ^6.5 GeV. 

To make comparisons of the measured dijet cross sections with the NLO QCD predictions, 
an asymmetric cut on the jet transverse energies of E^^^ > 8 GeV and E^^^ > 5 GeV 
was used for all further analysis. As seen in Fig. § E""^'""* = 8 GeV is sufficiently far 
away from the unphysical region near the symmetric-cut threshold. The predictions of 
NLO QCD with this cut, using DISENT with /i^ = Q^, agree with the data to within 
5%. This 5% difference between the predictions and the data is well within the estimated 
±15% systematic uncertainty of the measurements. 

The NLO QCD predictions using DISENT are also shown for the choice of = E^/4 and 
underestimate the measurements by 20%. Calculations were also done using MEPJET. 
The predictions of MEPJET agree with those of DISENT to within 5%. 



9.2 Measurements of the dijet cross sections 

The measured dijet cross sections for jets of hadrons as a function of log^Q Q^, logiQ^E^ / Q"^) 
and log^o ^ are given in Tables and are shown in Figs. In Fig. ||, the predictions 
of ARIADNE and LEPTO are compared to the data. The MC models are in general 
agreement with the shapes of the distributions but fail to describe the normalisation of 
the measured cross sections. Both models are LO and thus the uncertainties due to the 
choice of are substantial. 

The inclusive dijet differential cross section, da/dlog^Q Q^, for jets of hadrons is compared 
to the predictions of NLO QCD in Fig. |^(a). The measured differential cross section de- 
creases by two orders of magnitude in the range 

10 < Q2 < GeV^ and is well described 
by the predictions using = Q^. To examine deviations between the measurement and 
the predictions, the ratio of the measured cross sections and the predictions of DISENT 
using fij^ = and CTEQ4M is shown in Fig. ^(b). The NLO predictions are consistent 
with the measurements to within 10% for the entire range. This difference is attributed 
to the fact that the measured cross section includes events with two or more jets, whereas 
the NLO QCD calculation is limited to final states with two or three partons only[|. The 
study of inclusive dijets in this paper will permit comparison to higher-order calculations. 

The measured cross sections were also compared to the predictions of DISENT using 
= E^/4. For > 400 GeV^, no significant difference is seen between the predictions 
with = and with /i^ = E^/A. For lower values of Q^, the predictions with 
/x^ = E^/4: underestimate the measured cross sections by as much as 30%. Figure |^ 

^This hypothesis is supported by the fact that the exclusive dijet cross section (for two and only two 
jets) at high has been shown Q to agree to within 1% with the predictions of DISENT using = Q^. 
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also shows that the predictions of DISENT with = using the proton PDFs from 
MBFITIM [0 and those using the CTEQ4M PDFs agree equally well with the data. 
The NLO QCD calculations exhibit a significant uncertainty due to the choice of The 
scale uncertainty decreases from +20% at ~ 15 GeV^ to ±10% at = 400 GeV^ 
and ±5% at = 10^ GeV^. Only for Q'^ > 200 GeV^ are the theoretical uncertainties 
comparable to the systematic uncertainties of the data. 

To examine the interplay of the Q and Et scales, da / dlogif){E^/ Q'^) is compared to the 
predictions of DISENT in Fig. p. The predictions with yu|, = Q"^ using either CTEQ4M or 
MBFITIM agree well with the data. The scale uncertainty is ±(15 — 20)% in the region 
where < Q'^ and rises to ^20% for E^ > Q^. It is noteworthy that the prediction with 
= gives an excellent description of the data for all values of E^/Q^, whereas, as 
can be seen from Fig. |(b), the predictions using = E^/i fail to describe the data for 
Et > Q^. 



9.3 Gluon distribution of the proton 

Figure ^(a) shows da /dlogiQ ^. The requirement that two jets with high transverse energy 
be observed in the final state suppresses the cross section in the low-^ region. Therefore, 
the measured dijet cross section rises in the interval —2.3 < log^^Q^ < —1.7. For logj^Q^ > 
—1.7, the cross section decreases due to the decrease of the gluon and quark densities at 
high X. 

In Fig. ^(a), the measured cross sections are compared to the predictions of NLO QCD 
calculated with DISENT using both fij^ = and = E'^/A. The predictions with 
= Q2 are shown for both the CTEQ4M and the MBFITIM proton PDFs. Figure |(b) 
demonstrates that the NLO QCD predictions have large theoretical uncertainties over 
most of the ^ region. These uncertainties arise from the sensitivity to the choice of ^\ 
in the low-Q^ region, as shown in Fig. ^(b). Given these uncertainties, all of the pQCD 
calculations shown are consistent with the data. 

The measured dijet differential cross sections as a function of ^ for different ranges are 
given in Tables ^ and |^ and compared to the NLO QCD predictions in Figs. 10 and 11 



At low Q^, the uncertainty on the measurements is similar to the difference between the 
predictions using CTEQ4M and MBFITIM. As increases, the dijet cross section is 
limited by the kinematic and jet cuts to high ^ values. For > 215 GeV^, all of the 
predictions of NLO QCD give consistent results to ±5%, independent of the choice of /i^ 
iE^/4: or Q2) or proton PDFs (CTEQ4M or MBFITIM). The NLO renormalisation-scale 
uncertainty in this range decreases from 10% to 5% as increases. Over the entire 
range of these measurements, the predictions of DISENT with = are in reasonable 
agreement with the data. 

In the range 10 < < 464 GeV^, the total experimental uncertainty of the measure- 
ment ranges from (10 — 15)% and the theoretical uncertainty, due to the variation of 
varies from 45% to 10%. In the low-Q^ region, the dijet sample is dominated by 
gluon-induced BGF events and so could, in principle, be used to determine the gluon 
density. However, the size of the scale uncertainties precludes a precise extraction of the 
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gluon density^. The agreement between the data and the DISENT calculation indicates 
that the parametrisations of the gluon distribution in CTEQ4M and MBFITIM, which 
are determined primarily from the measurements of the proton structure function, are 
consistent with these dijet measurements. 

10 Conclusions 

Differential dijet cross sections in neutral current deep inelastic scattering have been 
measured in the ranges 10 < < 10^ GeV^, y > 0.04, E'^ > 10 GeV, > § GeV, 

E^f^ > 5 GeV, E^f^ > 5 GeV, EJ^f^ > 5 GeV and \r]^^^\ < 2 using the longitudinally 
invariant /ct cluster algorithm in the Breit frame. 

The measured cross sections were compared to next-to-leading-order QCD calculations 
as implemented in DISENT and MEPJET. The measured dijet cross sections fall by two 
orders of magnitude over the range considered here. The NLO QCD predictions agree 
with the data to within 10% when the renormalisation scale is chosen to be Q. The 
estimated theoretical uncertainties of up to 50% for < 20 GeV^ arise from the absence 
of higher-order terms, leading to a sensitivity on the choice of the renormalisation scale. 
Only for > 200 GeV^ are the theoretical uncertainties comparable to the systematic 
uncertainties of the data. The NLO predictions agree with the measured E^/Q"^ and 
^ distributions to within 10%. Again, for these distributions, the renormalisation-scale 
uncertainty is large; up to 40% when E^ > and 50% when log^^Q^ < —2. 

At low Q^, the uncertainty on the measurements is similar to the difference between the 
predictions using the CTEQ4M and MBFITIM parametrisations of the proton PDFs. In 
the range 10 < < 464 GeV^, the total experimental uncertainty of the measurement 
ranges from (10 — 15)% and the theoretical uncertainty, due to the variation of fij^, varies 
from 45% to 10%. Within these uncertainties, the measurements are consistent with a 
QCD prediction based on a gluon distribution extracted from scaling violations of the E2 
structure function of the proton. In this sense, therefore, these measurements are con- 
sistent with a universal gluon PDF. However, the large theoretical uncertainty precludes 
a useful determination of the gluon PDF in the proton using the dijet data; improved 
calculations are needed to exploit the full potential of this data. 
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logioQ^ bin 




Astat 


Agys 


Aes (pb) 


ISR/FSR 

correction 


1.00, 1.20 


1120.5 


± 18.5 


+154.3 
-141.4 


+92.3 
-93.0 


1.024 ±0.012 


1.20, 1.40 


967.7 


± 17.7 


+80.0 
-95.8 


+84.5 
-85.6 


1.040 ±0.014 


1.40, 1.60 


853.6 


± 16.7 


+46.0 
-66.9 


+68.5 
-75.0 


0.990 ±0.014 


1.60, 1.80 


690.5 


± 15.1 


+34.3 
-24.0 


+60.9 
-56.3 


1.002 ±0.016 


1.80, 2.00 


549.4 


± 13.3 


+35.0 
-14.6 


+45.5 
-44.2 


0.987 ±0.018 


2.00, 2.33 


368.2 


± 7.3 


+16.1 
-11.5 


+24.1 
-26.1 


0.989 ±0.012 


2.33, 2.67 


219.1 


± 5.1 


+5.3 
-5.3 


+13.8 
-14.0 


1.004 ±0.011 


2.67, 3.00 


113.5 


± 3.5 


+2.2 
-4.3 


+5.7 
-5.9 


1.021 ±0.012 


3.00, 3.50 


44.2 


± 1.6 


+2.4 
-0.6 


+2.2 
-2.4 


0.980 ±0.007 


3.50, 4.00 


8.7 


± 0.7 


+0.8 
-0.7 


+0.3 
-0.4 


0.960 ± 0.008 



Table 1: Dijet cross sections, dcr/rflog^o where is in GeV^, for jets of hadrons in 
the Breit frame, selected with the inclusive /ct algorithm in the ranges 10 < < 10^ 

GeV2, y > 0.04, E'^ > 10 GeV, E^^^ > 8 GeV, E^f^ > 5 GeV, E^^^ > 5 GeV, EJ^f^ > 5 
GeV and < 2. The statistical, systematic and jet- energy-scale uncertainties are 

shown separately. The multiplicative correction applied to account for effects of initial- 
and final-state QED radiation is shown in the last column. 
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Xog.^m/Q^) bin 



-1 


70, 


-1 


00 


-1 


00, 


-0 


70 


-0 


70, 


-0 


30 


-0 


30, 


+0 


00 


+0 


00, 


+0 


30 


+0 


30, 


+0 


70 


+0 


70, 


+1 


00 


+1 


00, 


+1 


70 



da/d\og,o{EyQ^) A, 



stat 



^sys 



Aes (pb) 



+1.8 


+0.9 


-1.2 


-0.6 


+3.9 


+4.3 


-2.4 


-6.7 


+7.5 


+9.4 


-5.2 


-11.4 


+21.6 


+27.9 


-12.6 


-26.6 


+30.7 


+46.3 


-26.7 


-48.8 


+45.0 


+75.7 


-45.6 


-73.5 


+14.6 


+64.7 


-30.4 


-67.4 


+1.9 


+11.9 


-6.4 


-11.9 



ISR/FSR 

correction 



27.6 
96.5 
213.2 
429.5 
640.9 
901.3 
688.7 
113.6 



± 1.1 
± 3.4 
± 5.0 
± 9.4 
± 12.2 
± 12.7 
± 12.5 
± 2.8 



0.943 ±0.010 
0.978 ±0.013 

1.009 ±0.012 
1.000 ±0.014 
0.994 ±0.013 

1.010 ±0.010 
1.021 ±0.013 
1.029 ±0.017 



Table 2: Dijet cross sections, da/dlogiQlE^/Q"^), for jets of hadrons in the Breit frame, 
selected with the inclusive algorithm. For details, see the caption of Table 1. 
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logio C bin 



-2.40, -2.20 
-2.20, -2.00 
-2.00, -1.80 
-1.80, -1.60 
-1.60, -1.40 
-1.40, -1.20 
-1.20, -1.00 
-1.00, -0.80 



dfj/dlogioC A 



stat 



»sys 



^ES 



198.2 
636.4 
1084.5 
1251.4 
1147.9 
653.7 
273.9 
82.6 



± 8.1 



-1-13.1 
-10.4 



±14.8 till 

±19.5 1:11-1 



± 20.6 



-1-63.3 
-63.1 



±20.0 till 



± 15.0 
± 9.5 
± 5.0 



-1-39.8 
-24.1 

+11.1 
-9.7 

-)-3.8 
-7.0 



-1-26.4 
-25.6 

-1-55.8 
-56.1 

-1-83.0 
-88.4 

+81.7 
-85.0 

+74.2 
-80.2 

+53.4 
-42.4 

+17.5 
-21.5 

+6.3 
-3.9 



ISR/FSR 



correction 



1.138 ±0.037 
1.052 ±0.018 
1.043 ±0.013 
0.995 ±0.011 
0.980 ±0.011 
0.994 ±0.013 
0.986 ±0.019 
0.997 ±0.031 



Table 3: Dijet cross sections, da/dlogi^^ ^, for jets of hadrons in the Breit frame, selected 
with the inclusive kr algorithm. For details, see the caption of Table 1. 
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logio C bin 


da/dlog 


LO C ^stat 


Asys A 


ES (pb) 


ISR/FSR 














correction 








1.0 < logio 


g2 < 1.2 






-2.40, 


-2.20 


70.2 


±4.9 


+7.0 
-3.9 


+6.6 
-7.2 


1.173 ±0.065 


-2.20, 


-2.00 


169.9 


±7.7 


+20.3 
-21.5 


+14.6 
-16.7 


1.021 ±0.034 


-2.00, 


-1.80 


261.9 


±9.9 


+25.9 
-27.8 


+21.2 
-19.9 


1.073 ±0.028 


-1.80, 


-1.60 


251.5 


±9.3 


+22.2 

-21.4 


+15.7 
-15.4 


1.002 ±0.024 


-1.60, 


-1.40 


199.6 


±8.4 


+ 18.8 
-21.7 


+13.9 
-16.9 


1.012 ±0.026 


-1.40, 


-1.20 


78.4 


±5.0 


+12.9 

-3.7 


+6.3 
-5.2 


0.967 ±0.034 








1.2 < logio 


< 1.4 






-2.40, 


-2.20 


55.9 


±4.6 


+1.5 
-7.1 


1 Q 1 

+8.1 
-8.8 


1.214 ±0.077 


-2.20, 


-2.00 


142.2 


±7.2 


+13.2 
-13.5 


+14.8 
-11.3 


1.092 ±0.041 


-2.00, 


-1.80 


222.5 


±9.3 


+7.3 

-13.3 


+17.8 

-19.6 


1.089 ±0.032 


-1.80, 


-1.60 


212.8 


±8.6 


+17.0 
-10.7 


+14.9 
-10.8 


1.003 ±0.027 


-1.60, 


-1.40 


181.2 


±8.2 


+9.4 
-14.9 


+ 11.2 
-15.1 


0.968 ±0.028 


-1.40, 


-1.20 


83.3 


±5.6 


+15.0 

-2.2 


+6.1 
-7.1 


1.108 ±0.042 








1.4 < logio 


< 1.6 






-2.40, 


-2.20 


33.9 


±3.2 


+3.8 
-2.4 


+5.9 
-3.2 


1.084 ±0.083 


-2.20, 


-2.00 


128.0 


±7.0 


+3.0 
-17.2 


+8.2 
-11.3 


1.064 ± 0.044 


-2.00, 


-1.80 


182.7 


±8.2 


+2.9 
-8.7 


+12.9 
-16.9 


0.988 ± 0.032 


-1.80, 


-1.60 


198.4 


±8.5 


+8.8 
-19.2 


+12.6 
-16.6 


0.938 ±0.028 


-1.60, 


-1.40 


165.2 


±7.9 


1 1 1 7 

-8.6 


J_1 ^ 

-r±o.O 

-9.6 


1.000 ± 0.031 


-1.40, 


-1.20 


89.5 


±6.1 


+5.1 
-6.7 


+7.9 
-9.5 


0.977 ± 0.040 








1.6 < logio 


< 1.8 






-2.40, 


-2.20 


25.5 


±2.9 


+1.2 
-1.4 


+3.2 
-3.9 


1.085 ±0.100 


-2.20, 


-2.00 


83.0 


±5.3 


1 1 7 
+ 1.7 

-7.7 


+8.7 
-7.6 


1.032 ± 0.050 


-2.00, 


-1.80 


146.4 


±7.3 


+11.6 
-8.8 


+11.5 
-12.1 


1.059 ±0.038 


-1.80, 


-1.60 


172.0 


±8.0 


+21.3 
-3.9 


+13.2 
-11.2 


1.000 ±0.033 


-1.60, 


-1.40 


149.1 


±7.6 


.1. Q Q 

-1.7 


_1_Q 7 

+y. / 
-10.8 


0.942 ± 0.032 


-1.40, 


-1.20 


69.2 


±5.2 


+3.5 
-6.0 


+8.4 
-3.9 


0.988 ± 0.046 








1.8 < logio 


< 2.0 






-2.20, 


-2.00 


61.6 


±4.5 


+2.6 
-5.7 


+6.5 
-4.7 


1.019 ±0.057 


-2.00, 


-1.80 


113.0 


±6.3 


+4.5 
-3.7 


+6.2 
-8.6 


1.008 ±0.041 


-1.80, 


-1.60 


133.7 


±6.9 


+8.7 
-3.3 


+12.2 
-9.6 


1.019 ±0.037 


-1.60, 


-1.40 


124.5 


±6.8 


+7.0 

-3.7 


+8.0 
-9.6 


0.941 ±0.035 


-1.40, 


-1.20 


66.8 


±5.1 


+7.2 
-5.2 


+8.1 
-5.2 


0.943 ± 0.044 



Table 4: Dijet cross sections, da/dlog^^^, for jets of hadrons in the Breit frame, selected 
with the inclusive kr algorithm. For details, see the caption of Table 1. 
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logio C bin 


I oiU 


^stat 


A sys 


Aes 


ISR/FSR 














mrrprtinn 






Z.U ^ lU 


glO V 


<^ Z.dO 






—2.20, 


—2.00 


43.2 


±3.2 


+5.3 
-1.7 


+2.3 
-3.2 


1 1 o« _|_ n nc,/i 
i.izu Hz u.uo^ 


nn 
— z.uu, 


1 8n 
— i.oU 


105.9 


±5.1 


+9.5 
-7.0 


+8.3 
-8.5 


n QQ7 4- n noQ 
u.yy ( It u.uzy 


— l.OU, 


1 fin 

— l.OU 


153.2 


±6.5 


+3.6 
-7.6 


+8.1 
-11.7 


980 ± 024 


— l.DL), 


1 Ad 

— 1.40 


148.7 


±6.5 


+6.4 
-3.1 


+9.3 
-7.1 


972 ± 023 


1 /I n 
— i.4U, 


1 on 
— i.zU 


94.7 


±5.4 


+2.2 
-6.6 


+6.3 
-5.8 


n Q81 -1- n no8 
u.yoi zt u.uzo 






2.33< logioQ^ 


9 fi7 
z.o / 






— z.UL), 


1 on 
—1.80 


44.2 


±3.0 


+2.3 
-3.3 


+4.6 
-2.3 


1 nqn _|_ n nqo 


— l.oL), 


—1.60 


87.1 


±4.5 


+3.0 
-1.9 


+2.5 
-5.6 


1 nfio -1- n no/i 


1 fin 

— l.DU, 


1 An 


104.0 


±5.1 


+3.9 
-5.2 


+5.2 
-7.6 


n Q88 -1- n non 
u.yoo zt u.uzu 


— 1.4(J, 


1 on 
— 1.20 


75.7 


±4.2 


+7.6 
-2.4 


+5.9 
-2.5 


n ofiA + n noi 


1 on 
— i.zU, 


1 nn 
— i.UU 


38.1 


±3.2 


+3.4 
-0.8 


+2.4 
-3.5 


n Q7Q -1- n nqn 
u.y / y zt u.uou 






2.67 < logio 


\ o.u 






— Z.UU, 


1 on 
— l.oO 


6.4 


±0.9 


+1.0 
-0.6 


+0.6 
-0.4 


1 nco _|_ n n«o 
i.uoz zn u.uoz 


— I.OL), 


1 i^n 

— I.dU 


31.1 


±2.2 


+0.4 
-1.2 


+1.4 
-2.5 


n QQo -1- n n9Q 
u.yyz zn u.uzy 


— l.DL), 


— 1.41) 


47.5 


±2.9 


+3.1 
-1.9 


+2.6 
-2.1 


1 ni n -1- n no/i 

I.UIU zt U.UZ4t 


— 1.4(J, 


— i.zL) 


54.2 


±3.3 


+1.7 
-4.9 


+1.6 
-2.1 


1 (YXA -I- n no/i 

i.Uo4 zt U.UZ'i 


— l.zL), 


—1.00 


33.1 


±2.6 


+2.7 
-2.1 


+2.6 
-1.9 


1 ncQ _|_ n noQ 
i.uoo zt u.uzy 






3.0<logiog' 


< 6.0 






1 on 
— l.oU, 


1 ^n 
— l.bO 


5.2 


±0.9 


+1.5 
-1.5 


+0.5 
-0.8 


1 ncn _L n nqo 
i.UDU zt U.Uoo 


1 /^n 
— l.bU, 


1 yl n 

— i.40 


18.7 


±1.7 


+2.9 
-1.0 


+0.6 
-1.4 


n Q81 -1- n ni 8 
u.yoi zt u.uio 


1 /in 
— i.4U, 


1 on 
— 1.20 


32.3 


±2.3 


+ 1.3 
-1.6 


+2.1 
-1.5 


n QQO -1- n ni 4 
u.yoz zt U.U1+ 


1 on 
— i.zU, 


1 nn 
— I.UU 


28.9 


±2.1 


+2.2 
-0.6 


+ 1.1 
-0.9 


n Q«o _|_ n ni /I 
U.yOo zt U.U1+ 


—1.00, 


— U.8U 


18.6 


±1.7 


+0.9 
-1.1 


+0.4 
-0.9 


n QfiQ -1- n ni 7 
u.yoy zt u.ui < 


n sn 

— U.oU, 


n fin 

— U.DU 


5.6 


±0.9 


+0.7 
-0.6 


+0.7 
-0.3 


n OS1 4- n nqo 
u.yoi zt u.uoz 






3.5<logioQ2 








-1.40, 


-1.20 


0.9 


±0.4 


+0.7 

-0.0 


+0.3 

-0.0 


1.092 ±0.041 


-1.20, 


-1.00 


6.4 


±1.0 


+0.6 
-1.1 


+0.0 
-0.5 


0.981 ±0.018 


-1.00, 


-0.80 


5.7 


±0.9 


+ 1.0 
-0.2 


+0.3 
-0.1 


0.957 ±0.014 


-0.80, 


-0.60 


6.9 


±1.0 


+0.3 
-0.7 


+0.1 
-0.3 


0.923 ±0.015 


-0.60, 


-0.40 


1.9 


±0.5 


+0.5 
-0.4 


+0.1 
-0.1 


0.964 ±0.027 



Table 5: Continuation of Table 4. For details, see the caption of Table 1. 
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Figure 1: (a) The QCD-Compton diagram and (b) the boson-gluon fusion diagram. The 
fraction of the proton momentum, p, carried by the stuck parton is ^. 
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Figure 2: Normalised uncorrected distributions of (a) logio(Q^/GeV^), (b) logio(-E'j'/Q^), 
(c) logio^ and (d) logiox for dijet events selected with the kr algorithm in the Breit 
frame, for E^^^ > 8 GeV, E^^^ > 5 GeV, E^^^ > 5 GeV, E^^^ > 5 GeV, \ri^^^\ < 2, 
y > 0.04 and > 10 GeV in the range 10 < < 10^ GeV^. The points are the 
measurements. The statistical uncertainties are generally smaller than the symbols. The 
solid and dashed histograms are the predictions of the ARIADNE 4.10 and LEPTO 6.5 
Monte Carlo programs, respectively. 
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Figure 3: Normalised uncorrected distributions of (a) E^^^, (b) E^^^, (c) ^^ab ^^^^ 
?72^^ for dijet events selected with the algorithm in the Breit frame. For details, see 
the caption of Fig. 2. 
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Figure 4: Contours representing a fixed scale uncertainty of the NLO QCD calculations 
as a function of log^o and logiQ(-E'|,/4), where Q"^ and E'j. are in GeV^. The scale 
uncertainty shown in (a) was calculated by taking the ratio of the predicted dijet cross 
sections cr(/x|j = / 4) / a {nj^ = Q^); that shown in (b) was calculated by taking the ratio 
of the predicted dijet cross sections cr(/i^ = E'^/IQ) / a{ii\ = E^/A). The cross sections 
were calculated using DISENT and refer to jets of partons selected with the inclusive kx 
jet algorithm in the Breit frame. CTEQ4M was used for the proton PDFs. 
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Figure 5: The dijet cross section for jets of hadrons in the Breit frame, selected with the 
inclusive kx algorithm as a function of E'^^^''^"*. The points represent the measured cross 
section. The typical magnitude of the statistical and systematic uncertainties added in 
quadrature is ±15%. The full line represents the predictions of NLO QCD using DISENT 
with = Q^. The dashed line is the prediction of DISENT with = The dotted 

line is the prediction of MEPJET with = All three pQCD calculations use the 

CTEQ4M proton PDFs, refer to jets of hadrons and were obtained by dividing the parton- 
level predictions of DISENT by correction factors computed using the ARIADNE 4.10 
MC simulation. 
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Figure 6: Dijet cross sections (a) da/d\ogiQ{Q'^/GeV'^), (b) da/d\ogiQ{E^/Q'^) and (c) 
da /dlogiQ ^ for jets of hadrons in the Breit frame selected with the inclusive fcy algorithm. 
The points represent the measured cross sections. The error bars are generally smaller 
than the points. The shaded band represents the systematic uncertainty due to the jet- 
energy scale. The full (dashed) histogram represents the predictions of ARIADNE 4.10 
(LEPTO 6.5) with the CTEQ4M proton PDFs. 
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Figure 7: (a) Dijet cross section, da/dlogiQ Q^, for jets of hadrons in the Breit frame se- 
lected with the inclusive kr algorithm. The points represent the measured cross sections. 
The inner error bars represent the statistical uncertainties and the outer bars are the 
statistical and the systematic uncertainties added in quadrature. The shaded band repre- 
sents the systematic uncertainty due to the jet-energy scale. The full line represents the 
predictions of NLO QCD using DISENT with = and the CTEQ4M proton PDFs. 
The dashed-dotted hne is the prediction of DISENT with = E^/A and the CTEQ4M 
PDFs. The dashed line is the prediction of DISENT with = using the MBFITIM 
proton PDFs. (b) The cross sections in (a) divided by the predictions of DISENT with 
A*ii ~ CTEQ4M PDFs. The cross-hatched band represents the theoretical 

uncertainty due to the choice of calculated by choosing = Q'^/A and /^^ = Q'^. 
All of the pQCD predictions shown here refer to jets of hadrons and were obtained by 
dividing the parton-level predictions of DISENT by the average value of the correction 
factors computed using the LEPTO 6.5 and the ARIADNE 4.10 MC simulations. The 
shaded band in (c) shows the magnitude and the uncertainty of this parton-to-hadron 
correction. 
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Figure 8: Dijet cross section, da/dlogiQ{E^/Q'^), for jets of hadrons in the Breit frame 
selected with the inclusive algorithm. Other details are as described in the caption of 
Fig. 0. 
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Figure 9: Dijet cross section, da/dlogiQ^, for jets of hadrons in the Breit frame selected 
with the inclusive kj- algorithm. Other details are as described in the caption of Fig. 0. 
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Figure 10: Dijet cross sections, da/dlogi^^ ^, for jets of hadrons selected with the inclusive 
kx jet algorithm in different regions of Q^. The points represent the measured cross 
sections. The inner error bars represent the statistical uncertainties and the outer bars 
are the statistical and the uncorrelated systematic uncertainties added in quadrature. 
The shaded band represents the systematic uncertainty due to the jet-energy scale. The 
full line represents the predictions of NLO QCD using DISENT with = and the 
CTEQ4M proton PDFs. The dot-dashed line is the prediction of DISENT with = 
and the CTEQ4M proton PDFs. The dashed hue is the prediction of DISENT with 
= g2 ^gjj^g ^Yie MBFITIM proton PDFs. The cross-hatched band at the bottom of 
each plot is the /i^ uncertainty calculated by varying //^ in the range (4Q^, (5^/4). All of 
the pQCD predictions shown here refer to jets of hadrons and were obtained by dividing 
the parton level predictions of DISENT by the average value of the correction factors 
computed using the LEPTO 6.5 and the ARIADNE 4.10 MC simulations. 
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Figure 11: Continuation of Fig. for higher values of Q^. Details are as described in the 



caption of Fig. 10. 
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